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ABSTRACT 


A theoretical  analysis  is  presented  of  unsteady  solid 
propellant  combustion,  particularly  combustion  stability  .and 
extinction  by  rapid  depressurization.  It  is  assume^*  that 
the'  solid  decomposes  by  a pyrolosis  law  and  the  gasdous  prod* 
ucts  react  exothermically  following'  ah  Arrhenius  law.  For 
large  values  of  this  non-dimensional  activation  energy  the 
gas-phase  combustion,  turns  out  to  be  quasiateady.  The-  char 
acteristic  response  time  of  the  solid  to  gas-phase  perturb- 
ations turns  out  to  be  large  compared  to  the  characteristic 
residence  time  in  the  heat-up  zone  of  the  solid,  their  ratio 
being  of  the  order  of  the  non-dimensional  activation  energy 
in  the  gas-phase.  A linear  stability  criterion  has  been  ob- 
tained that  gives  stable  burning  for  steady  burning  rates  a- 
bove  and  below  two  limiting  valuesjthe  width  of  the  unstable 
region'  increases  with  the  activation  energy  of  the  pyrolisis 
law,  and  becomes  zero  for  a finite  value  of  that  activation 
energy. nonlinear  stability  analysis  has  been  performed 
that  sl/ows  the  existence  of  a limit  cycle  for  unsteady  burn- 
ing conditions j in  a part  of  the  cycle  the  burning  is  slow 
and  the  soVid  responses  as  a whole  to  variations  in  gas  phase 
burning;  duVing  other  part  the  characteristic  time  is  short 
and  the  heat\content  of  the  solid  is  constant.  The  dynamic 
extinction  process,  has  also  been  investigated. 
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N ,0  M E N C L A..T  U R E 


A Amplitude  of  limit  cycle,  eq.(64), 

B Preexponential  factor,  and  constant  of  integration  (eq.52). 

c Specific  heat. 

C Constant  of  integration  (eq.52). 

D Diffusion  coeficient.  j 

I 

E Activation  energy  of  the  pyrolisis.  i 

\ ' J 

II 

E Activation  energy  of  the  gas-phase  reactioVx.  ! 

. O 

E'  Dimensionless  activation  energy  of  the  gas-phase  reaction. 

1 Dimensionless  heat  flux  at  the  surface,  see  Eq..  (^14). 

L Heat  of  vaporization. 

L Lewis  number, 

e 

m-  Burning  rate. 

n Exponent  of  fuel  mass  fraction  in  chemical  reaction  rate, 

p Pressure. 

q Total  heat  qonteiit  of  the  solid  (eq.  71). 


P. 


% 
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Total  heat  content  of  the  solid  (eq.67). 

Q Heat  released  per  unit  mass  of  fueli 
R Universal  gas  constant. 

T Temperature. 

t Time. 

t Characteristic  time,  see  Eq.(43). 
c 

u Velocity 
X Space  coordinate. 

Y Fuel  mass  fraction. 


y Dimensionless  space  coordinate. 

a Thermal  difussivity,  and  dimensionless  parameter  ('eq.62). 

Y Dimensionless  steady  state  surface  temperature,  see  Eq.(38). 


6 Dimensionless  width  of  the  unstable  region.  Eq.(57). 


e 

e 

X 


Dimensionless  inverse  activation  energy 
Nondimensional  temperature. 

Thermal  conductivity. 


B n 
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y’  Dimensionless  burning  rate. 


. y'  Dimensionless  burning  rate  (eq.  58). 


i Nondimensional  space  coordinate. 


IT  Npndimensional  pressure  history. 


p Density. 


Dimensionless  time,  see  Eq.(19). 


T Dimensionless  time,  see  Eq.(30). 


T Dimensionless  time,  see  Eq. (61). 


(I)  See  Eq.(38)  and  (46). 
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f Flame. 


g Gas. 


s Surface. 


o Steady  state  conditions. 


-»  Infinity  in  the  solid.. 
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1 . Introduction  ^ 

An  understanding  of  the  unsteady  burning  of  solid  propel- 
lants provides  insight  into  such  important  problems  as  the  pres 
surization,  depressurization  and  stability  of  a solid-rocket 
motors.  ^ 

Typical  mo4els  used  for  the  analysis  are  those  presented 
in  references  (1)  to  ((12), 

Research  in  the  stability  of  combustion  in  a rocket  motor-, 
has  been  active  over  the  past  twenty  years,  since  irregular 
pulses  in  chamber  pressure  were  observed  to  develop  instead 
of  the  expected  smooth  pressure-time  history,  these  irregular 
pulses  are  generally  accompanied  by  more  regular,  small-ampli- 
tude, pressure  oscillations,  with  frecuencies  of  the  order  of 
the  natural  vibrational  frecuencies  of  sound  waves  in  the'  cham- 
ber. Combustion  instability  leads  to  inefficient  operation  of 
rocket  motors  and  even  to  mechanical  failure  of  the  propellant. 

A considerable  theoretical  effort  has  been  devoted  towards 
the  understanding  of  this  phenomenon.  However,  there  are  , so  many 
different  effects  that  may  influence  the  stability  of  a -burning 
solid  as  to  prevent  a consensus  on  the  theoretical  description 
of  the  phenomenon. 

An  excellent  review  of  the  analyses  of  the  small  amplitude 
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pressure  oscillatjions  has  been  developed  by  Culick  . These 
analyses  calculate  the  admittance  function  of  the  burning  surface 
which  gives  the  burning  rate  response  to-  a pressure  disturbance, 
indicating  therefore  whether  pressure  oscillations  are  am;  lified 


Qr  attenuated.  Most  of  the  analyses'  assume  the  gas  phase  to  be 
quas'isteady , in  the  sense  that  the  gas  phase  adjusts  very  quick- 
ly to  changes  in  conditions  when  compared  to  the  response  of  the 

^ ' 

solid,  and  differ  mainly  in  the  assumptions  used  to  calculate 
the  heat  transfer  from  the  gas  to  the  solid  phase.  However, as 
Culick  points  out,  the  majority  of  the  results  lead  to  the 
same  two-parameter  form  of  the  admittance  function,  with  different 
definitions  of  the  two  parameters.  A stability  boundary  in  the 
space  which  coordinates  the  two  parameters  entering  the  admittance 
function,  is  defined  as  the  curve  where  the  admittance  function 
becomes  infinite,  so  that  a small  pressure  change  causes  a large 
fluctuation  in  burning  rate.  Points  above  this  stability  boundary 
produce  unstable  solutions  and  points  below  this  boundary  provide 
stable  steady  solutions. 

Interest  in  the  response  of  a burning  solid  to  an  externally 

applied  pressure  variation  stems  from  the  possibility  of  extinction 

3 

by  a rapid  pressure  decay  . This  dynamic  extinction  process  is 
useful  in  order  to  design  solid  propellant  rockets  with  ctcp- 
restart  capabilities,  of  interest  in  connection  with  space  appli- 
cations . 

Nearly  all  the  theoretical  models  developed  to  explain  dy- 
9-14 

naraic  extinction  invoke  some  kind  of  quasisteady  ap- 
proximation for  the  gas-phase.  However,  as  pointed  out  in  ref- 
erence  (12)  some  of  these  analyses  have  interpreted  incor- 

rectly this  assumption  by  implying  that  the  heat  feebback  from 
the  flame  to  the  solid  is  a steady  state  function  of  the  instan- 
taneous pressure  only. 


The  present  paper  is  an  attempt  towards  analyzing  unsteady 
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processes  in  solid  propellant  burning  by  means  of  asymptotic 
techniques  based-  on  the  assumption  that  the  nondimensional 
activation  energy  of  the  gas-phase  reaction  is  large.  Attend 
tion  i*s  paid  both  to  the  problem  of  the  stability  of  the 
steady  state  and  to  the  response  of  the  burning  propellant  to 
an  externally  imposed  pressure  variation. 

We  consider  a one-dimensional  model  in  which  a condensed 

material  gasifies  by  a rate-controlled  surface  process  and  then 

reacts  in  the  gas  phase.  This  gas  phase  reaction  is  described 

by  an  Arrhenius  law,  and  we  consider  the  limit  in  which  the 

nondimensional  activation  energy  is  large.  The  quasisteady 

assumption  is  used  to.  describe  the  gas-phase,-  so  that  we  may 

3 is 

use  the  results  obtained  by  Williams  and  Buckroaster  et  al  ,. 
when  analyzing  the  quasisteady  burning  of  a solid  in  the  limit 
of  high  activation  energy.  These  analyses  yield  the  burning 
rate  and  the  heat  feedback  to  the  solid  as  functions  of  the 
pressure  and  the  flame  temperature,  and  these  relations  are 
then  used  to  analyze  the  unsteady  response  of  the  condensed 
phase.  It  is  found  that  the  characteristic  response  time  of 
the  solid  is  large, of  the  order  of  the  nondimensional  activation 
energy,  so  that  the  temperature  profiles  are  quasisteady  in 
first  approximation  and  the  transient  term  is  a perturbation  of 
the  quasisteady  solution  which  may  be  calculated  from  the  qua- 
sisteady profiles.  In  this  way  a differential  equation  is 
derived  which  describes  the  evolution  with  time  of  the  burning 
rate  as  a function  of  the  nondimensional  surface  temperature, 
the  nondimensional  pyrolisis  activation  energy  and  the  pres- 
sure-time history.  The  analysis  of  this  equation  yields  the 
stability  condition  of  the  steady  state  solution.  It  is  found 

that  burning  conditions  are  stable  outside  a range  of  values  of 
steady  burning  rates.  The  amplitude  of  this  unstable  range 


increases  with  the  activation  energy  of  the  pyrolisis,  and  becomes 
zero  for  a finite  value  of  that  activation  energy.  It  turns  out 
that  instability  is  related  to  the  fact  that  the  total  heat  con- 
tent of  the  solid  is  a decreasing  function  of  the  burning  rate 
under  steady  conditions. 

Then,  the  existence  of  non-linear  oscillations  for  condi- 
tio*\s  corresponding  to  instability  under  steady  conditions  is 
investigated.  However,,  the  equation  obtained  for  the  burning 
rate  shows  that,  as  the  unstability  boundary  is  crossed,  the 
time  derivative  becomes  infinity  revealing  that  during  part  of 
the  oscillation  period  there  is  a characteristic  time  much  shorter 
than  the  one  corresponding  to  unsteady  response  of  the  solid.  The 
previously  obtained  equation  is  conveniently  modified,  and  then 
it  is  found  that  the  nonlinear  oscillations  due  to  instability 
are  amplified  until  they  reach  a limit  cycle  that  is  character- 
ized by  two  intervals  of  time;  one  of  then  is  long  and  during 
it  the  temperature  profiles  in  the  solid  are  quasisteady,  in 
the  other  the  characteristic  time  is  much  shorter  and^  the  solids 
has  no  time  to  change  its  heat  content.  The  amplitude  of  this 
limit’  cycle  increases  with  the  activation  energy  of  the  pyrcHisls. 
For  Values  of  this  activation  energy  close  to  the  one  that 
corresponds  to  the  disappearance  of  the  unstability  region,  the 
amplitude  of  the  limit  cycle  becomes  zero,  and  the  period  of 
this  cycle  decreases,  disappearing  the  difference  between  the 
two  time  intervals,  now  the  variation  of  the  burning  rate  with 
time  is  smooth  and  in  the  limit  tends  to  be  sinusoidal. 

The  dfrferential  equation  describing  the  quasisteady 

evolution  of  the  burning  rate  with  time  is  then  used  to  obtain 

w 

the  response  of  the  solid  to  a pressure  decrease  and  a possible 
criterion  for  dynamic  extinction  is  obtained. 
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2.  Formulation 

We  consider  a one-dimensional  model  with-  the  solid  occu- 
pying the  half  space  x < 0 and  the  gaseous  phase  the  region 
X > 0.  For  convenience,  the  origin  is  fixed  at  the  surface  of 
the  regressing  solid.  Figure  1 is  a schematic  representation  of 
the  process  and  shows  the  effects  accounted  for  in  this  study. 

Surface  gasification  is  assumed  to  occur  by  an  irrever- 
sible pyrolisis  process  which  is  described  by  an  Arrhenius  law. 

A one  step  over-all  exothermic  reaction  takes  place  in  the 
(premixed)  gas.  The  present  model  has  been  used  because  it  has 
been  successful!  in  describing  the  steady-state  deflagration  of 
several  propellants. 

The  equation  of  conservation  of  momentum  reduces  to  the 
•statement  that  the  pressure  is  approximately  uniform  throughout 
the  region  treated  but  varies  with  time^,  and  we  will  alsjq  use 
the  well  justified  assumptions  that  the  work  assoc, iated  with 
viscous  and  external  forces  is  negligible  and  we  will  use  Fick’s 
law  to  calculate  the  diffusion  velocities. 

With  those  assumptions,  the  conservation  equations  of  mass, 
concentration  and  energy  in  the  condensed  phase  and  in  the  gas 
phase,  become  respectively. 
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where  all  symbols  are  defined  in  the  nomenclature. 


These  equations,  together  with  the  equation  of  state  should 

determine  p , u,  Y,  T,  T . 

. ® , S 

The  boundary  conditions  are 


T(t,  o'*')  = T(t,  O")  = T 


. 3T  3T 

^ 3x  s"  3x 


|gt  " + L 
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T - T at  x-»'-«,Y=0,T  = T£  atx->-»  (8) 

n * 

The  pyrolisis  process  is  assumed  to  follow  an  Arrhenius  law 


m = B’  exp(-B/RT  ) 

s 


where  m is  the  mass  flux  relative  to  the  burning  surface,  and  E 
the  activation  energy  of  the  pyrolisis  process. 

3.  Gas-phase  analysis 


In  most  studies  of  unsteady  solid  propellant  burning,  the 

gas-phase  is  assumed  to  be  quasisteady  in  the  sense  that  the  ]jes-  ' 

ponse  time  in  the  gas  is  short  compared  to  the  response  time  of 

the  condensed  phase.  The  ratio  of  these  characteristic  times  is 

of  the  order  of  the  ratio  of  the  thermal  responsivities  of  solid 

2 

and  gas,  which  is  usually  small  i The  quasisteady  assumption  will 
not  be  valid  when  analyzing  the  response  of  a burning  solid  to 


h 

.4 
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Very  high  frecuency  pressure  oscillations  . 


When  the  gas-phase  is  considered  to  be  quasisteady,  the  ana- 
lysis is  greatly  simplified  since  the  conservation  equations  are 
uncoupled.  All  the  time  derivate  terms  may  be  neglected-  so  that 
equation  (1)  reduces  to  p^u  = m.  In  addition  the  dp/dt  term  in 
equation  (3)  may  be  neglected.  Under  these  conditions  the  gas-phase 
equations  may  be  solved  in  the  limit  of  high  activation  energy  of  the 
gas-phase  reaction.  This  quasisteady  solution  has  been  derived  by 

3 

Williams  in  the  case  of  Lewis  number  unity  and  by  Buckmaster  et  ■ 
al^^ j for  arbitrary  Lewis  number,  Le . 

Rather  than  repeat  those  analyses  we  will  only  state  the  result^ 
needed  in  the  following  sections.  The  reader  is  referred  to  reference 

I' 

3 and  15  for  details.  The  analysis  in  reference  (15)  considers  Le  cons- 
tant but  arbitrary.  In  the  limit  of  high  activation  energy,  it  is 
found  that, the  reaction  term  is  only  significant  in  a thin  flame 
where  T is  close  to  T^^.  In  terms  of  the  nondlmensional  variables 


8 = 


5 = / (m/a)dx 


the  solution  outside  of  the  flame  sheet  becomes  simply 


- 1 + 1 exp(0 

8 

1 - exp{L^(5  - Cj)} 


for  5 < 


(11) 


for  K > Kf  " (j 

y = 0 

where  the  location  of  the  flame  sheet,  is  given  to  leading 

order  by 

= In  Cl/1)  (: 


M 


The  parameter 


1 = I 
^ dC  s 


is  the  dimensionless  heat  conducted  out  of  the  gas  at  the  interface. 
The  effect  of  the  solid  and  the  pyrolisis  are  concealed  in  1*  An 
overall  energy  balance  in  the  gas  phase,  proyi.des  a relationship 
between  0 , 0^  and  1. 

Si 


1 t 


(15) 


Since  the  flame  sheet  must  lie  in  the  gas , equation  (13) 
provides  limitations  in  the  possible  values  of  1 


0 < 1 r 1 + 0g  - 


(16) 


When  1 0 the  flame  sheet  moves  to  infinity  and  no  heat  reaches 

the  condensed  phase  from  the  flame.  When  1 1 the  flame  sheet 

approaches  the  surface  and  all  the  heat  generated  at  the  flame  goes 
to  the  solid. 

Solutions  for  1 -»■  0 and  1 1 are  presented  in  reference  (15) 

The  reaction  zone  is  located  in  the  vicinity  of  where 
(5  - is  of  order  0|^/E’^.  The  parameter  E'  is  the  nondimensional 

activation  energy  E'  = E c /RQ,  In  the  reaction  zone  the  teltnpera.-: 
ture  differs  from  the  flame  temperature  by  a small  quantity  of  order 
T-R/E  . To  leading  order,  only  the  reactive  and  diffusive  terms 
are  important  in  this  zone.  Solution  of  the  energy  equation  with 
the  appropriate  matching  conditions  to  the  frozen  solution  outside 
of  the  flame  sheet,  provides  the  burning  rate  eigen~value. 

In  this  way  an  expression  for  the  burning  rate  is  obtained 


14.  - 


Involving  the  pressure  and  the  flame  temperature,  namely 


m = 2r(n+l)aB0^’^^^/E'"'^^  exp(- 


For  large  activation  energy  the  effect  of  the.  exponential 
term  is  dominant  so  that  the  square  root  term  may  be  taken  as 
constant  when  analyzing  small  changes  in  the  flame  temperature. 

5 

Equation  (17)  coincides  with>  the  Denison-Baum  formula  , 


4.  Condensed  Phase  Analysis 

The  characteristic  time  in  the  condensed  phase  is  short 

2 

compared  with  the  characteristic  time  in  the  gas  phase  . Therefore, 
the,  condensed  phase  should  not  be  considered  quasisteady.  However, 
a brief  description  of  the  quasisteady  solution  will  be  presented 
before  considering  the  unsteady  analysis. 

In  terms  of  the  nondimensional  variables  for  the  soil'd 


T -T 
so  «» 


m c 

y = 


■ 2 
m_c 

^1  ""  W ^ 


M = ^ 
m 

o 


Eq.  (‘4)  may  be  written 


ay  . ay 


(20) 


The  subscript  o refers  to  the  initial, 
steady  condition.  The  boundary  condition  (6)  may  be  written, 


:i5.- 


I 


by  using  Eqs.  (14)  and  (15)  as 


o_e 

9y 


c(T„^-T  ) 

so  oj ' 


cT„  - + Q - L + ce^(T^^-T  ) 

" g r . s so  <» 


(21) 


To  solve  the  quasisteady  problem,  we  may  neglect  the  time 
derivative  term  in  Eq.  (20),  The  solution  to  Eq.  (20)  is  simply 

• ■ 6=6^  exp  (yy)  (22) 

S 

so  that  from  Eq,  (21). 

^g^fo  = + Q - L (23) 

% 

This  equation  can  also  be  obtained  by  an  overall  energy 
balance  and  shows  that  the  flame  temperature  is  constant  in  a 
quasisteady  process.  For  a given  T^,  Eq.  (23)  gives  the  value  of 
T^^  which  can  be  used  in  the  Ibiison-Baum  formula,  Eq.  (17),.  to 
calculate  m.  The  pyrolisis  law,  Eq.  (9), yields  T , so  that  Eqs. 
(11),  (12)  and  (22)  describe  the  complete  temperature  profile. 

Figure  (2)  shows  the  burning  rate  m under  steady  conditions; 
as  a function  of  pressure  for  a fixed  value  of  T^  and  therefore 
of  Tj^.  This  curve  is  calculated  by  using  Eqs.  (17)  and  (23). 

The  surface  temperature  T is  a parameter  along  the  curves 

S 

of  figure  (2),  since  it  is  related  to  tho  burning  rate  m^  through 

the  pyrolisis  law.  However,  T^  is  limited  by  the  inequalities  (16) 

so  that  only  a portion  of  the  curve  applies.  At  P„  the  flame 

sheet  has  moved  off  to  infinity.  At  P the  flame  sheet  has  reached 

s * 

the  surface. 


i 


If  wmiwiM 
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We-  will  now  analyze  the-  evolution-  with:  time  of  an  initially 
steady  temperature  profile  during  an  unsteady  process.  After  some 
manipulations,  Eqs . (20)  and  (21)  provide  the  following  relations- 
hips 


i®. 

3y 


s 


(24) 


where  use  has  been  made  of  Eq.  (23).  The  left  hand  side  of  the 
preceding  equation  represents  the  difference  betv/een  the  heat 
flux  existing  during  an  unsteady  process  and  the  one  existing  if 
the  process  was  quasisteady.  ' We  will  consider  unsteady  processes 
that  result  in  changes  of  order  unity  in  the  burning  rate  with 
respect  to  the  one  existing  under  quasisteady  conditions.  Equa- 
tion (17)  shows  that  in  the  limit  of  high  activation  energy  of 

the  gas  phase  reaction,  small  changes  in  the  flame  temperature  of 
2 

order  RTi/E  , produce  variations  of  order  unity  in  the  burning 

* O 

rate.  Therefore,  the  right  hand  side  of  Eq.(24)  is  small,  so 
that  the  heat  conducted  to  the  solid  during  an  unsteady  process 
differs  by  a small  quantity  from  the  heat  conducted  if  the  process 
■was  quasisteady. 


The  second  equality  of  Eq.(24)  indicates  that  the  variation 

with  time  of  the  heat  content  in  the  heat  up  zone  of  the  solid 

2 

is  small,  of  order  RT^/E  . When  this  heat  content  decreases,  the 
solid  appears  from  the  gas  phase  as  a heat  source,  and  therefore 
the  flame  temperature  increases. 


It  i 
changes  of 
the  solid, 
therefore 


s necessary  to  wait  times  of  order  E^/RT^  to  produce 
order  unity  in  the  heat  content  of  the  heat  up  zone  in 
The  characteristic  response  time  of  the  solid  is 
long  compared  with  the  characteristic  residence 


17.- 


time  in  the  heat  up  zone. 


Th^!  Denison-Ba.um  formula,  E^.  (1.7),  yields 


■ ^fo  ' *’o  - 2RT, 


which  in  the  limit  of  high  values  of  E /2RT^  becomes 


(26) 


g fo 


where 


" ,"g<V^fo>, 

- = exp  r } 

ir  2RT^^ 


it= 


Let’s  define  a small  parameter  e,  as 


c 2rt;: 

e - ■ --  fo 

c E (T  -T  ) 
g so  «' 


and  introduce  as  nondimensional  time  variable 

T = ex^  (30) 

Equation  (27)  may  be  used  to  express  Eq.  (20)  and  the 
boundary  condition  (24)  in  terms  of  the  new  time  variable  x,  as 


30  ^ 30 

' 37  37 


- e In 
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The  value  of  the  surface  temperature  may  be  written  in 
terms  of  the  burning  rate  y,  through  the  pyrolisis  law,  Eq.(9), 
resulting : 


_ Y-1 


lny+1 


s 1-wlny 


where 


Introducing  the  expansions 


8=6-+  eO . + ... 
0 1 


y = y-  + ey,  + 


Q T t-Hi  T . 


(35) 

(35'') 


in  Eqs.  (31)  and  (32)  the  following  equations  defining  6^  are 
obtained 


3y 

3y^ 

> 

3y 

= >'o®os 


(36) 


Whose  solution  is  the  quasisteady  solution 


% = ®so 


Integrating  eq.(31)  with  boundary  condition  (32)  the  following 
integral  condition  is  obtained; 


If  dy  = -y  in  ^ 


if 


19?" 


On  using  the  first  approximation  for  Sj  the  integral  in 
the  left  hand  side  may  be  calculated,  thus  obtaining 


d 

dT 


sO 

y 


-Poln 


It 


(39) 


where  Qgo/^^O  ^sP’^^sents  the  total  heat  content  of  the  solid 
under  steady  state  conditions.  Now  the  pyrolisis  law  is  used 
to  relate  the  sui’face  temperature  to  the  burning  rate,  obtaining 

<iVio  (y-1)  (l~&)lny-)y^  y 

= ^ - r ■■  In  (-2)  (HO) 

dT  Yw-K’+l+u(Y“2)lnyQ+w  In  Pq  it 


which  describes  the  evolution  with  time  of  the  burning  rate 

during  an  unsteady  process,  in  terms  of  two  parameters  y 

to.  For  a quasisteady  evolution  dy^/dx  vanishes  and  therefore 

Pq  = TT. 


The  value  of  the  burning  rate,  p^,  obtained  from  Eq.(HO), 
may  be  used  in  the  pyrolisis  law  to  deduce  the  surface  temper- 
ature history,  and  through  Eq.(37)  the  complete  temperatui>e 
profile . 

5 . Stability  Analysis. 

In  this  section  an  analysis  is  presented  of  the  stability 
of  the  steady  deflagration  of  a solid  which  undergoes  an 
Arrhenius  type  gas  phase  reaction  with  large  activation  energy. 

• « 

The  pressure  is  considered  to  be  constant,  so  that  ir  = l. 

The  burning  rate  of  the  steady  solution  is  m^ , so  that  initially 
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y-=l.  Let's  assume  that  at  t=0  a perturbation  changes  the  value 
of  the  burning  rate  so  that  yQ  = l+y'.  Eci.(40)  shows  that 


iH.’-  _ 

dT 


Therefore  the  stability  of  the  steady  state  solution  depends 
on  the  value  of  the  parameter  • If  this  parameter  is  greater 

than  one,  the  solution  is  stable,  being  unstable  when  it  is  smaller 
than  one.  In  the  stable  case  the  perturbations  decrease  exponen- 
tially. The  parameter 


E(T  ^-T  ) 


is  identical  v/ith  the  parameter  A used  by  Denison  and  Baum'^  when 
analyzing  the  unstable  burning  of  solid  propellants.  These 
authors  develop  a linearized  analysis  of  the  response  of  a burning 
solid  to  a pressure  disturbance  using  a model  which  parallels  the 
one  used  in  the  present  analysis.  The  difference  however,  is  that 
in  our  analysis  we  retain  the  nonlinear  effects  that  were  line- 
arized in  reference  (5).  In  figure  (1)  of  their  paper,  Denison- 
Baum  show  the  stability  boundary  which  separates  regions  of 
stable  from  regions  of  unstable  burning.  Two  parameters  A and  a, 
define  the  stability  of  a given  solution.  The  parameter  a is 
essentially  identical  to  our  parameter  e,  and  in  that  figure  (1) 
we  observe  that  for  a 0 the  solutions  are  stable  for  A < 1 and 
unstable  for  A > 1.  Therefore  our  stability  criterion  coincides 
with  the  one  of  Denison  and  Baum. 


The  instability  of  the  one-dimensional  deflagration  may  be 


clearly  interpreted,  when  the  activation  energy  of  the  pyrblisis 
is  large  (to  small).  Let's  assume  that  at  a certain  time  a per- 
turbation causes  the  burning  rate  to  increase  (decrease)  with 
respect  to  its  steady  value.  Since  the  activation  energy  of  the 
pyrolisis  is  large,  the  surface  temperature  will  remain  nearly 
constant.  However,  the  width  of  the  heat  up  zone  in  the  solid: 
will  decrease  (increase),  so  that  the  total  heat  content  of  the 
solid  will  decrease  (increase).  The  variation  of  the  thermal 
energy  of  the  solid  produces  an  increase  (decrease)  in  the  flame 
temperature,  as  may  be  seen  from  Eq.  (24),  which  will  further 
increase  (/decrease)  the  burning  rate.  This  selfaccelerating 
behaviour  results  in  instability  of  the  one-dimensional  defla- 
gration. When  the  activation  energy  of  the  pyrolisis  is  not 
large,  there  are  two  effects  which  control  the  response  of  the 
burning  solid  to  a perturbation  in  the  burning  rate.  When  the 
burning  i^ate  increases  the  width  of  the  heat  up  zone  decreases 
and  the  surface  temperature  increases.  The  decrease  in  width 
tends  to  decrease  the  heat  content  of  the  solid  thus  producing 
instability,  and  the  increase  in  surface  temperature  tends  to 
increase  the  heat  content  of  the  solid  thus  producing  stability. 
The  stability  of  the  solution  depends  on  which  of  these  two 
effects  dominates. 

If  the  reference  pressure  is,  such  that  ir  is  differ- 
ent from  one,  the  perturbation  of  the  steady  solution  will  be 

= TT  + y' 


and  Eq.(40)  will  show  that 


^ap8 


22 


where  F(Tr)  is  the  first  factor  on  the  left  hand  side  of  Eq.(4p) 
with  TT  = Hq  . An  alternative  form  of  writing  (43)  is,  from’ 
Eq.(39) 


■ll-O 


dy  * 


(43* ) 


When  the  stability  of  other  points  besides  the  one  that  corre- 
sponds to  ir  = l is  studied,  it  is  necessary  to  analyze  the  sign 
of  the  first  factor  on  the  right  hand  side  of  Eq(43’)‘;  when  it 
is  positive  the  steady  state  is  stable.  As  it  may  be  deduced 
from  Eqs.(38)  and  (43'),  the  steady  state  is  stable  when  the 
total  heat  content  of  the  solid  (under  steady  conditions) 
increases  as  the  burning  rate  is  increased;  this  criterion  is 
valid  for  any  pyrolisis  law. 

For  the  pyrolisis  law  given  by  Eq.(33)  the  stability  is 
related  to  the  sign  of  F(tt)  in  Eq.('43)  or  the  first  factor  of 
the  left  hand  side  of  Eq.(40).  In  this  factor  the  numerator 
vanishes  for  the  limiting  values  of 

y^  = exp(l/a))  and  Wq  = 0 (44) 


which  correspond  to  0 tending  to  infinity  (eq.(33)),  and 

s 

the  burning  rate  tending  to  zero  respectively.  It  will  be  ' , 
considered  that  y^  is  outside  these  limits.  The  stability 
criterion  is  then  given  by  the  sign  of  the  denominator  that 
vanishes  for: 


lnyo=  - 


^ Y-2  + /yCy-1/  w) 


(45) 


tuesacti 


lii'itfte'ilfiyfinrgHiiftVinf^ 


23.- 


For  values  of  ^0  = 7^  above  the  larger  root  or  below  the  lower 
root  of  Eq.(45)  the  steady  state  solution  is -stable. 

There  is  a certain  arbitrariness  in  the  choice  of  the 
value  of  P^. , (Eq.  (28)  ),  that  establishes  the  steady  state 
solution  for  Tr  = l . For  the  following  analysis  it  is  more 

Then  from  Eq.(38). 

(46), 

Equation  (45)  becomes 

InWo  = i ^ /iTau  (47) 

With  this  choice  of  the  reference  value  of  T^^ , the  steady 
state  solution  that  corresponds  to  11  = 1 is  unstable  and  is 
located  in  the  middle  of  the  unstable  region  (if  we  consider 
Iny  as  the  appropiate  variable). 

For  values  of  the  activation  energy  of  the  pyrolisis. 
smaller  than  a finite  value: 

E-<  4RT  , u > 1/2 

OO 

there  are  no  unstable  states.  The  thickness  of  the  unstable 
region  increases  as  the  activation  energy  of  the  pyrolisis 
increases. 

« 

The  pyrolisis  law,  Eqs.(33)  and  t9),  in  the  new  variables 


convenient,  to  choose  T =2T  . 

so  » 


2RT. 


Y = 2 


becomes : 
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6.  Non  linear*  stability  analysis 


The  purpose  of  this  section  is  to  show  that  unstable 
steady  states  produce  non-linear  oscillations  of  finite  amplitude 
around  the  steady  state.  Equation  (40),  however,  is  not  suf- 
ficient to  describe  these  oscillations.  A small  perturbation 
around  the  steady  state  will  start  to  grow  exponentially  as 
shows  Eq.(43)  (F<0).  As  this  perturbation  becomes  larger,  the 
linear  analysis  fails  and  equation  (40)  should  be  used.  This 
Eq.  (40)  gives  a monotonic.increase  of  y^,  and  when  it  gets  close 
to  the  value  given  by  Eq.(47)  its  derivative  becomes  infinite 
and  the  present  analysis  fails. 


Let  us  . etain  in  Eqs.(35  and  35')  the  terms  of  order  e 
and  it  will  be  shown  that  these  terms. will  introduce  a second 
order  derivative  of  y in  equation  (40),  which  generates  oscilla* 
tions  of  y with  time.  Additionaly,  to  retain  the  terms  of  order 
e will  make  the  first  derivatives  of  y finite  in  its  whole 
range  of  variation,  and  it  will  be  shown  that  Eq.(40)  will  be 
valid  only  for  some  intervals  of  those  oscillations. 

Retaining  terms  of  order  e in  Eq.(35)  and  ubing  Eq.(37). 

e=0g  exp(yy)  + eO^  (49) 


where  6^  should  satisfy  the  conditions 


6^  = 0 for  y = 0 and  y 


25.  T 


» 

■Introducing  Eq.(49)  into  Eq.(31)  and  neglecting  terms  of  order 
2 

E , the  following  equation  is  obtained 

de  dn  30.  3^0. 

e exp(vy)(^  + + ev  ^ =e  (51) 

The  solution  of  Bqs.  (50)  and  (51)  is; 


= (By  + Cy^)  exp  yy  (52) 

where 

1 d0  0 dy 

B = - — ^ — • ■ (53) 

y dT  y d'c 


(54) 


‘ This  solution  should  also  satisfy  the  boundary  condition 
(32).  As  it  was  done  previously,  this  condition  is  substituted 
by  an  integral  condition  obtained  by  integrating  Eq.(31)  between 
- » and  0,  and  using  Eq.(32). 

i \ ff"  + O(e^)  (55) 

•"CO 

Where  the  last  term  is  included  because  in  the  gas-phase 

. 2 

analysis  terms  of  order  e were  neglected.  By  introducing 
Eqs.(49),  (52),  (53)  and  (54)  into  (55)  the  following  equation 
for  y is  obtained  : 
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Where  8 is  a function  of  y given  by  the  pyrotisis  law  (eqs.(37) 
or  (48)),  although  in  principle  Eq. (56)  is  valid  for  any  pyroii- 
sis  law.  If  terms  of  order  e are  neglected,  this  .Eq.  becomes 
the  same  as  Eq.(39)  or  Eq.(40).  In  principle  the  term  of  order  e 
on  the  right  hand  side  which  has  n^t  been  calculated,  and  is  due 
to  higher  order  contributions  of  the  gas  phase  reaction,  is.  as 
important  as  the  new  calculated  term  (the  first  one),  unless, 
as  it  is  in  the  present  case,  the  time  scale  is  very  small  in 
regions  where  the  first  term  should  be  retained.  ‘The  last  term 
on  the  right  hand'  side  represents  the  influence  of  higher  o'rder 
terms  that  would  appear  if  this  procedure  would  be  continued. 

6.1.  Case  for  which  the  instability  region  is  narow 

Before  considering  the  more  general  case,  let  us  study 
the  case  for  which  « is  close  to  0.5. 

• « 

« = I (1-6)  (57) 

' » 

6 <<  1 


Then,  from  equation  (47)  vhe  instability  range  will  be 


Iny  = + 2/6 


Since  we  are  going  to  study  values  of  y located  near  ox'  in  the 
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instability  range, we  make  the  change  of  variables 


y = 1 + 2y ' 


(58) 


The  point  of  steady  burning  corresponds  to  = ir  and 

.we  assume  this  point  to  be  in  or  near  the  instability  region 


If  = 1 + 2y'/? 

o 


(59) 


Then,  equatibn  (56)  in  the  new  variables  becomes,  after  using 
equation  ^(48) 

\ 


-e6 


1/2 


d^y » 
dT^ 


26®^^'.(y’^“l)^ 

dt 


26^/^(;y’' 


(60) 


+ 0(  ed 


d^y’ 

dT^ 


) + 


0(6^  iEl)  + 

dt 


0(e)  . jt^)  = 

dT^ 


In  this  case  the  significant  limit  corresponds  to  a choice 
of  the  new  time  scale  that  will  be  shorter  than  the  previous 
one 


T = T /2/e  . - ( 61 ) 

and  6 should  be  small,  of  order  /e  • Introducing  the  parameter  a 


6 = d /e  / 2 


(62) 
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T,hen,  equation  (60)  becomes: 


dx^  dx 


dy ' 


(y’  ^ y;)  = 


(63) 


and  all  remaining  terms  become  negligible.  Equation  (63)  is  a 

Lienard  type  equation,  that  reduces  to  Van  der  Pol's  equation  when 

y'=0.  For  |y'|>l  the  steady  state  solution  is  stable-  since  we 
o 'o' 

are  outside  the  instability  region.  For  ly^|<l  the  steady  state 
is  unstable  for  a small  perturbation;  however,  as  the  perturba- 
tion grows  it  tends  towards  a limit  cycle  of  finite  amplitude. 

This  means  that  if  we  maintain  a constant  pressure  that  would 
correspond  to  a steady  state  burning  in  the  unstable  region,  t'he 
burning  rate  would  exhibit  an  oscillating  behavior  of  finite 
amplitude. 

If  y^  is  close  to  the  instability  boundary,  (y^  close  to 
one)  the  amplitude  of  the  limit  cycle  will  be  small  and  the  so- 
lution of  (63)  for  the  limit  cycle  will  be  of  small  amplitude, 
and  can  be  aproximated  by: 


y ' = y ' + A sen  x 
o 


(64) 


Then  to  impose  the  condition  that  the  secular  ti^rms  of  , 
equation  (63)  will  vanish,  let  us  multiply  eq.(63)  by  - 

and  integrate  through  a cycle: 


o 


6V\2 
d-  X ^ 


dT  = 


0 
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pbtaining 

A = (65) 

Results  (64)  and  (65)  for  the  limit  cycle  are  valid  only  fp,r 
a small  or  j | close  to  one. 

When  a is  large,  the  solution  approaches  that  which 
corresponds  to  an  instability  region  of  finite  amplitude, 
since  from  equation  (62')  6 will,  become  larger  (see  also  eq.(57)). 
Now  one  may  explore  the  limit  a large  in  eq. (63),  (still 
maintaining  the  limit  6 small). 

— 

Rewritting  eq.(63)  in  the  form  (defining  a new  variable  q)^ 
a ^ =y’-y’  (66) 


1 ^ = q 

o dx 


(67) 


and  analyzing  the  phase  plane  corresponding  to  equations  (6'6) 
and  (67)  for  and  a large  (see  fig. 3a) 


dy 


y;-y’ 


q-(H 


I 3 


y') 


(68) 


q is  the  total  heat  content  of  the  solid.  The  meaning 
of  this  new  variable  and  the  physical  interpretation  of  the 
results  are  reserved  for  the  next  section,  where  they  will 
become  more  clear. 


'“*r 
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it  can  be  seen  that  if  the  steady  burning  point  is  in  the 
unstable  branch  BD,  then  it  will  go  out  from  that  point  and. 
approach  the  branch  AB  or  CD  and  from  there  it  would  finally 
approach  the  limit  cycle  ABCD.  This  limit  cycle  has  an  ampli^ 
tude  (corresponding  to  points  C and  A)  of: 

y’  = ,i  2 

While  we  are  in  branches  AB  or  CD  in  the  stable  region,  the 
solution  is  given  by: 

a 4:  - V’)  = y'-y'  (69) 

'dT  3 ° 


which  is  esentially  equation  (40)  corresponding  to  the  case  in-- 
which  the  second  derivative  of  y with  respect  to  time  is  negli- 
gible. Then,  the  characteristic  time  is  long;  tT  is  of  order  a. 
In  branches  BC  and  DA  on  the  other  hand  "q'-'  remains  approxi- 
mately constant,  equal  to  t 2/3,  then,  the  characteristic  time 
is  short  of  order*  T'vl/a,  for  the^  variation  of  y,  given  now  by 
equation  (67).  ' 

In  figure  3b  the  phase  plane  is  presented  for  the  case 

in  which  y'  is  in  the  stable  region, 
o 

This  method  for  a largo  feiibs  when  l.y^  |.  is  sufficiently 
close  to  !■;  then  the  limit  cycle  would  be  given  by  equations 
(64)  and  (65).  Figure  *4  presents  the  results  of  the  numerical 
integration  of  equation  (63);  giving  the  maximum  and  minimuin 
Value  of  y'  in  the  limit  cycle  for  different  values  of  y^  and  a 
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6 . 2 Case  for  which  the  instability  region  is.  finite. 

To  extend  this  analysis  to  the  case  when  the  instability 
rang^e  -has  a thickness  of  order  unity,  the  method  presented  on 
sectim  6.1  for  o large  will  be  used.  Assuming  equation  (56) 
to  be  valid,,  and  rewritting  it  in  a form"  .similar  to  that  of 
equations  (66)  and  (67),  the  following  relations  are  obtained 


^ = 
dT 


in 


V; 

7 


, d(0  /y^)  dy 
•t  s ^ s 

y dy  dx  y 


(70) 


(71) 


From  equations  (55)  and  (70)  it  can  be  seen  that  "q"  represents 
. the  total  heat  content  of  tho  solid.  As  long  as  the  time  x is 
of  order  unity,  the  first  terra  of  the  right  hand  side  is  small 
and,  q=0_/y,  as  it  happened  for  the  quasisteady  solution 
(equation  (39)).  However,  if  during  the  limit  cycle  there  is 
a period  in  which  x've  (that  is  Xj^''*!,  the  characteristic  time  . 
for  heat  conduction  in  the  solid,  see  equations-  (19)  and  (30)), 
then  the  total  heat  content  of  the  solid,  q,  would  he  given 
by  the  complete  equation  (71).  Nevertheless,  equation  (.55) 
expresses  that  the  characteristic  time  to  transfer  across  the 
solid  surface  the  heat  needed  to  change  the  total  heat  content 
of  the  solid  is  of  order  x.  In  order  words  if  equation  (70) 

is  rewritten  with  the  characteristic  time  x,  it  becomes: 

1 
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that  means  that  q remains  almost  constant  during  periods  of 
order  during  this  period  there  is  no,  time  for  transfer  of 
the  heat  needed  to  change  the  total  heat  content  of  the  solid. 

Figure  (5)  shows  schematically  0 /y  as  a function  of'y, 

s 

and  by  analogy  with  figure  3a  one  may  assume  that  the  limit 
cycle  is  given  by  ABCD.  In  branch  AB  the  process  is  quasi- 
steady and  corresponds  to  making  e equal  to  zero  in  equation 
(71),  or  using  Eqs.(39)  or  (40).  In  branch  BC  q is  constant, 
and  the  characteristic  time  is  t^,  y changes  very  rapidly  :and 
there  is  no  time  for  the  heat  content  of  the  solid  to  change., 
this;  process  continues  until  we  are  well  in  the  stable  region 
in  the  point  where  y has  the  value  that  makes  under  steady 
conditions  Figui'e  5b  shows  schematically  the  way  in 

which  y will  change  as  a function  of  time. 

The  previous  qualitative  results  are  independent  of  the 
pyrolisis  law.  We  are  now  going  to  calculate  the  maximum  and 
minimum  values  of  y (v^»  Pq)  tor  the  pyrolisis  law  expressed 
by  equation  (48). 

For  point  A of  figure  5 ; 

. i„„)  = 1 (72) 

(l  + v^l-2u)  (1-ulny ) u 

For  point  C it  can  be  shown  that 
iny^  = -Iny^ 
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The  solution  of  eq.(72)  depends  only  on  parameter  u and 
is  presented  in  figure  6. 


For  u close  to  0.5 


ulnjJ^  = -2/i-2ti) 


that  corresponds  to  the  solution  of  section  6.1  for  a large. 
The  other  limiting  value  of  w corresponds  to  a large  value  of 
the  pyrolisis  activation  energy,  then  w tends  to  zero,  and 


Iny^-*-  -1/u; 


winy,  = -1  + -r- y exp(- 

^ (1-/1 -2  ur 


This  means  that  for  u small 


lnp->-il  /to 


and  from  equation  (48); 


e -»■  0,  «» 

s 


Obviously,  then,  the  pyrolisis  law  is  probably  not  correct, 
however,  this  is  an  indication  that  the  solid  is  igniting  and 
extinguishing  alternatively  in  its  oscillations. 


In  Eq.(72)  is  probably  not  well  justified  to  neglect 
terms  such  as  0(e2  4^)  of  Eq.(56),  that  will  originate  terms 
of  order  g in  the  right  hand  side  of  eq(71)»  This  is 
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ju^tified  in  the  analysis  of  section  6.1  when  the  unstability 
region  is  thin.  Then,  the  fast  transition  regions  corresponding 
to  branches  BC  and  DA  of  figure  5 will  not  be  appropiately  de- 
scribed by  Eq.(72)  with  q constant.  A different  approach  will 
be  needed  that  will  probably  consist  in  calculating  more  accu- 
rately the  total  heat  content  of  the  solid  when  v and  0^  change 
rapidly  (subject  to  the  pyrolisis  law,  eq.(48))  and  impose  the 
condition  that  this  heat  content  does  not  change  in  branches 
BC  and  DA.  • 


7 . Response  of  a burning  .solid  to  pressure  variations 


In  a previous  rspprt,  ref  (16),  corresponding  to  this 
grant,  an  analysis  was  presented  of  the  response  of  a solid 
to  an  imposed  pressure  variation;  the  evolution  with  time 
of  the  burning  rate  was  calculated  by  using  Eq.(40)  with 
the  non-dimensional  pressure,  n,  being  a known  function  of 
time  defined  by  Eq.(28).  Here  we  shall  only  reproduce  the 
main  conclusions  obtained  in  that  report,  ref.  ( 16  ). 

It  was  shown  in  section  4 that  the  characteristic  response 
time  of  the  solid  is 


Xp  E (T  -T  . ) 
■ ^ g so  ” 

' " "‘^f 


(74) 


Therefore,  if  the,  characteristic  time  of  pressure  variation 
is  small  compared  with  it,  the  pressure  variation  is  a step 
function.  If  the  time  of  pressure  variation  is  long,  compared 
to  the  characteristic  time  shown  in  Eq.(74),  the  response  of 
"the  solid  may  be  considered  as  a sequence  of  stationary  states. 
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^ In  the  present  analysis  only  depressurization  problems 
will  be  cphsidered,  so  that  n will  decrease  from  some  initial 
value  TTj  at  time  0,  to  its  final  value  iTp.  In  addition  it  is 
assumed  that  the  initial  solution  is  stable: 

either  where  = exp(^  /l-2w)  (75) 

or  < Hg  where  = exp(“  ) (76) 

where  and  1:2  are  the  values  of  the  pressure  that  limit  the 
unstable  region,  Eq.(47)-. 

Two  possible  types  of  evolutions  are  possible  a)  and  b):: 

a)  The  unstable  region  is  not  crossed"^  during  depressurizatipvi 
This  is  the  case  when: 

al)  The  unstable  region  does  not  exist 

4RT  > E ’ ■ 

V> 


a2)  The  initial  and  final  values  of  pressure  are  both  above 
the  unstable  region. 


a3 ) Both  values  of  pressure  are  below  the  unstable  region 


^^2  " 


^I  ^ 
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In  this  case  a)  the  evoUtion  of  the  burning  rate,  y, 
given  by  Eq.(40)  is  smooth  and  it  tends  monotonically  to  its 
final  steady  value  ifp  = Pp. 

b)  The  unstable  region  is  crossed  during  depressurization. 

. This  is  the  case  when 

bl)  Hj.  > > iTp  the  final  steady  state  is  stable. 

b2,)  TTp  > Tt^  > ifp  > Tfg  the  final  steady  state  is  unstable. 

In  this  case  b)  Eq.(40)  has  a singularity  at  y = ahd  at  a 
time  such  that  ttCt)  •-  ir^.  The  character  of  the  singularity 
depends  on  the  value  of  the  shape  of  the  pressure  time  function 
curve  when  the  instability  region  is  entered  at  When 

1 exp(~/l-2u) 

< (77) 

2u  /l-2u 


dir 

dt  1 


this  singularity  has  a node-like  character  and  Eq.(40)  is  able 
to  give  a transition  across  the  singular  p.oint.  On  the  other 
hand  when  condition  (77)  is  not  satisfied >the  singularity 'has 
a spiral  character  and  the  solution  ceases  to  be  valid  when  y 
approaches  since  the  variation  of  burning  rate  with  time 

is  so  rapid,  that  is  no  longer  valid,  to  assume  that  the  deriv- 
ative term  is  only  a perturbation  of  rhe  quasistoady  solution., 
A different  type  of  analysis , (may  be  similar  to  that  presented, 
in  section  6 of  this  report),  will  be  required  to  study  the 
behaviour  of  the  solution  near  y = if^.  However,  since  y is 
decreasing  very  rapidly  to  zero,  the  condition  oposite  to  (77)= 
-may  be  identified  as  a criterion,  for  dynamic  extinctipn  caused 


\ 


ntummr' 


37  i - 


by  a rapid  depressurization  process. 

When  condition  bl  is  satisfied  there  is  another  singular 
point  at  Vi  = Tf2  is  always  a saddle  through  which  the  solu- 

tion wiri  exit  from  the  unstable  region. 

Any  perturbation  occurring  while  the  value  of  y is  in 
the  unstable  region  will  grow  with  time  leading  to  the  non- 
linear oscillations  studied  in  section  6.  It  is  possible  that 
this  study  on  depresurizatipn  is  not  appropiately  described  by 
Eq.(40),  and  should  be  remade  taking  into  account  higher  order 
unsteady  effects  similar  to  those  retained  in  section  6 when 
studying  non-linear  oscillations. 

8 . Conclusions 

An  asymptotic  analysis  has  been  developed  of  the  unsteady 

burning,  of  a sold  propellant  in  the  limit  of  high  activation 

energy  of  the  gas  phase  reaction.  The  analysis  shows  that  the 

ratio  of  the  characteristic  response  time  of  the  solid,  and 

the  characteristic  residence  time  in  the  heat-up  zone  is  of  the 

ord^e.r  of  the  nondimensional  activation  energy  of  the  gas-phase 

reaction.  This  slow  response  of  the  solid  allows  to  suppose 

that  the  temperature  distribution  is  in  a first  approximation 

/ 

quasisteady,  and  by  integrating  the  condensed  phase  energy 
equation,  an  equation  is  derived  whose  solution  yields  the 
evolution  with  time  pf  the  burning  fate. 

It  is  found  that  for  values  of  the  activation  energy  of 
the  pyrolisis  sufficiently  large; 

E > 4RT 

00 
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there  are,  for  other  parameters  fixed,  a range  of  non-dimensional 
pressures,  eq.(28),  such  that 

- - A -2(0  < lniT<  - /l-2w  , 

U 44 

,in  it, steady  burning  will  be  unstable.  This  criterion  corresponds 

5 

to  that  of  Denison-  y Baum  in  the  limit  of  high  activation  energy 
of  the  gas  phase  reaction. 

The  evolution  with  time  of  the  burning  rate  for  a fixed 
pressure  corresponding  to  instability  has  been  investigated; 
the  non-linear  equation  obtained  from  the  quasisteady  temperature 
distribution  is  not  appropiate  to  study  this  evolution  because 
it  gives  derivatives  of  the  burning  rate  with  respect  to  time 
that  are  infinite  when  the  unstable  region  is  crossed.  This 
suggests  that  during  part  of  the  evolution  of  the  burning  rate 
the  appropiate  time  scale  should  be  shorter.  Higher  order 
terras  are  then  retained  when  calculating  the  temperature  distri- 
bution in  the  solid  and  a new  equation  for  the  burning  rate 
evolution  is  obtained;  this  new  equation  contains  second  order 
derivatives  of  the  burning  rate  which  g nerate  oscillations  of 
the  burning  rate  with  time.  This  new  equation  is  first  solved 
in  the  limit  that  the  instability  interval  is  small;  and  it  is 
found  that  when  the  width  of  this  interval  is  of  the  order  of 
fourth  root  of  the  inverse  tof  the  non-dimensional  activation 
energy  of  the  gas-phase,  thei-c  is  a significative  limit  in  which 
all  the  retained  terms  are  of  thp  same  order  and  the  expansion- 
procedure  is  justified.  Then,  instability  leads  to  oscillations 
that  correspond  to  a limit  cycle  that  is  described  by  an  equation 
of  the  type  of  Van  der  Pool  equation  slightly  modified.  When 
the  width  of  the  instability  region  is  much  smaller  than  the 
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above-mentioned  limit,  the  oscillations  have  a sinusoidal 
character,  and  when  it  is  much  larger,  two  different  time 
intervals  appear  during  the  oscillations  period,  one  in- 
terval is-  large  and  the  other  short.,  corresponding  this 
last  one  to  a sudden  jump  in  the  burning  rate.  Based  on  this 
last  result , the  limit  cycle  is  investigated  for  the  case  ' 
that  width  of  the  instability  region  is  finite  and  the  cor- 
responding oscillations  are  of  the  form  shown  schematically 
in  fig.  5b.  During  part  of  the  oscillation  period  the  char- 
acteristic time  is  of  the  order  of  the  characteristic  re-' 
sponse  time  of  the  solid  and  much  larger  than  the  char- 
acteristic time  for  heat  transfer,  the  temperature  profile 
in  the  solid  is  quasisteady,,  and  the  burning  rate  changes 
slowly  until  it  reaches  the  instability  regions  then,  the 
characteristic  time  becomes  much  shorter,  of  the  order  of 
the  characteristic  time  for  heat  transfer,  and  the  temper- 
ature distribution  in  the  solid  is  no  longer  quasisteady, 
Howeverj during  this  short  time, which  is  much  shorter  than 
the  characteristic  response  time  of  the  solid> there  is  no 
time  for  changing. the  total  heat  content  of  the  solid.  In 
the  intervals  BC  and  DA  of  fig.  5a  the  burning  rate  changes 
very  rapidly,  maintining  the  total  heat  content  of  the  solid, 
and  goes  across  the  unsteady  region  until  it  is  found  in  the 
oposite  stable  region  a burning  rate  that  would  correspond  to 
a quasisteady  temperature  distribution  in  the  solid  of  the 
same  total  heat  content. 

The  expansion  procedure  indicated  in  this  work  is  hot  ^ 
appropiate  to  describe  the  structure  of  the  fast  transition 
regions  BC  and  DA  of  fig.  5a  (for  a finite  width  of  the 
instability  region),  because  in  this  regions  the  problem 
in  the  soiid  is  esentially  unsteady,  and  any  method  that  would 
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consider  the  time  derivatives  of  the  temperature  as  a ^perturbation 
in  the  energy  equation  of  the  solid,  would  not  be  valid.  The 
appropiate  method  would  consist  in  solving  this  energy  equation 
subject  to  the  condition  that  the  total  heat  content  of  the 
solid  is  constant,  when  the  burning  rate  is  changing. 

The  analysis  presented  in  this  work  makes  use  of  a pyrolisis 
law  of  the  form  given  by  eq.(9)(or  eq.(48),  in  non-dimensional 
variables).  However,  the  method  only,  uses  this  law  for  specific 
applications  and  is  in  principle  valid  for  any  ether  law  that 
would  relate  the  burning  rate  with  the  surface  temperature. 
Instability,  as  shown  in  fig.  5a  would  correspond  to  an  interval 
of  burning  rates  in  which  the  total  heat  content  of  the  solid, 
under  steady  conditions,  would  decrease  as  the  burning  rate 
increases. 

There  are  other  possible  mechanisms,  besides  the  one 
examined  here,  that  would  prevent  the  growth  of  the  oscillations 
such  as  three  dimensional  effects. 

The  analysis  has  also  been  applied  to  study  the  response 
of  a burning  solid  subjected  to  a depressurization  process, 
and  assuming  that  the  steady  solution  at  the  initial  pressure 
is  stable.  For  this  analysis  it  was  used  Eq.(40)  that  corre- 
sponds to  quasisteady  temperature  distribution  in  the  solid. 

It  is  found  that  for  depressurization  processes  such  that 
the  initial,  intermediate,  and  final  pressure,  would  all  cor- 
respond to  steady  stable  burning,  the  variation  of  burning  rate 
with  time  is  smooth  and  monotonous,  even  foil  step-like  changes  . 
in  pressure,  although  in  this  last  case  the  burning  rate  does 
not  adjust  instantaneously  to  the  new  pressure,  and  approaches 


it  in  a long  time.  The  absence  of  dynamic  extinction  ap.pears 
to  coincide  with  the  suggestion  of  T'ien  ■ , who  considers  that 
heat  losses  is  the  mechanism  responsible  for  dynamic  extinction* 
Nevertheless,  our  analysis  seems  to  indicate  that  there  is  a 
condition  that  may  be  identified  with  the  onset  of  dynamic 
extinction  if  during  the  process  of  depressurization  the  unstable 
region  is  crossed.  A mathematical  singularity  occurs  when,  during 
the  depressurization  process,  the  instability  region  .is  entered j 
and  the  character  of  this  singularity  depends  on  the  slope  of 
imposed  pressure-time  curve  at  the  singularity.  When  this  sdope. 
is  larger  than  a certain  critical  value  defined  in  Eq.  (77),  the 
singularity  has  a spiral-like  character  that  originates  deriv- 
atives of  the  burning  rate  with  respect  to  time  that  are  infinite. 
This,  and  the  fact  that  any  small  perturbation  will  be  amplified 
in  the  instability  region,  indicate  the  need  of  retaining  in 
some  way  (which  may  be  similar  to  the  one  used  for  studying 
non-linear  oscillations)  unsteady  effects,  in  the  temperature 
distribution  within  the  solid  when,  during  the  depressurization 
process,  the  instability  region  is  crossed. 


42.- 


' . references 

1.  Williams,  F.A.,  ’’Response  of  a Burning  Solid  to  Small-Ampl atude  - 
• Pressure  Oscillations”,  J.  Appl,  Phys . , 33,  3153-3166  (1962). 

2.  Summerfield,  M. , Caveny,  L.H.,  Battista,  R.A.,  Kubota,  N. , 
Gostintsev,  Yu.  A.  and  Isoda,  H.  , ’’Theory  of  Dynamic  Extinguish- 
ment of  Solid  Propellants  with,  Special  Reference  to  Nonsteady 
Heat  Feedback  Law”,  AIAA  paper  70-667  (1970).. 

3.  Williams,  F.A.  ’’Quasi-Steady  Gas-Phase  Flame  Theory  in  Unsteady 
Burning  of  a Homogeneous  Solid  Propellant”  j,  AIAA  J.,  11, 
1328-1330  (1973). 

4.  Hart,  R.W.  and  Mc.Clure,  F.T.,  ’’Combustion  Instabi’laty : Acoustic 

Interaction  with'  a Burning  Propellant  Surface”,  J.  Chem.  Phys. 
Vol.  30,  p. 1503(1959).  • , _ 

5.  Denison,  M.R.  and  Baum,  E.,  ”A  simplified  Model  of  Unstable 
Burning  in  Solid  Propellants”,  ARS  Journal,  31,  1112-1122 
(1961). 

6.  Green,  L.  Jr.,  Jet  Propulsion,  26,  655  (1956). 

7.  Culick,  F.E.C.,  ”A  Review  of  Calculations  for  Unsteady  Burning 
of  a Solid  Propellant”,  AIAA  J.,  6,  2241-2255  (1968). 

8.  Ciepluch,  C.G.  j,  ’’Effect  of  Rapid  Pressure  Decay  on  Solid"''  . 
Propellant  Combustion”,  ARS  Journal,  31,  1584-1586  (1961). 

9.  Von  Elbe,  G.,  ’’Theory  of  Solid  Propellant  Ignition  and  Response 
to  Pressure  Transients”,  Interagency  Solid  Propulsion  Meeting 
3,  95  (1963). 


- in.  .■.M...'mnT»ii«iiii«ww>i^^  .1 .; ' ^ 

43.- 

i 

10^  Paul,  B.E.,  Lovine,  R.L.  and  Fong,  L.Y.,  "A  Ballistic  Expiana- 
tion  of  the  Ignition  Pressure  Peak",  AIAA  paper  64-121  (1964). 

11.  Horton,  M.D.,  Bruno,  P.S.  and  Graesser,  E.C.,  "Depressurization 
Induced  Extinction  of  Burning  Solid  Propellant",  AIAA  J.,  6 
292-297  (1968). 

12.  Merkle,  C.L.,  Turk,  S.L.  and  Summerf ield',  M,  , "Extinguishment 
of  Solid  Propellants  by  Depressurization:  Effects  of  Propellant 
Parameters",  AIAA  paper  69-176  (1969). 

13.  T’ien,  J.S.,  "A  Theoretical  Criterion  for  Dynamic  Extinction 
of  Solid  Propellants  by  Fast-Depressurization'*,  Combi  Sci.  6 
Tech.,  9,  37-39  (1974). 

14.  Baliga,  B.R.  and  T'ien,  J.S.,  "A  Theory  of  Flammability  Limits 
of  Solid  Propellants  at  Low  Pressure",  AIAA  paper  75-229  (1975). 

15l  Buckmaster,  J.D.,  Kapila,  A.K.,  and  Ludford,  G.S.S.  "Linear 
Condensate.  Deflagration  for  Large  Activation  Energy".  Acta 
Astronautica , 3,  593-614  (1976). 

16.  Crespo,  A.  and  Kindelan  M.  "Unsteady  Processes  in  Solid 

Propellant  Combustion"  Scientific  Report  1.  AFOSR  76-3049. 

/ 

(May  1977). 


L . 1ST  OF  FIGURES  _ 

Figure  1 Schematic  representation  of  the  combustion  process. 

Figure  2 Steady  state  burning  rate  as._.a'  function  of  pressjur.e. 

Figure  3a  Phase  plane  showing  instability  and  the  form  in  wliich 
the  limit  cycle  is  reached  when;  th'S  instability  re- 
gion is  thin,  a large,  and  the  steady  state  is  uns-  ^ 
table.  Eqs(66)  and  (67). 

Idem,  when  the  steady  state  is  stable. 

Amplitude  of  the  limit  cycle  when  the  instability 
region  is  thin,  for  different  values  of  the  Iteady 
burning  rate,  and  different  values  of  a.  Eq.(63)i. 

Total  heat  content  of  the  solid  under  steady  state 
conditions  as  a function  6f  the  burning  rate.,  and 
limit  cycle. 

Schematic  showing  the  evolution  of  burning  rate 
with  time,  during  the  limit  cycle,  when  the  width 
of  the  instability  region  is  finite. 
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Amplitude  of  the  limit  cycle  when  the  width  of  the 
instability  region  is  finite'.  Eq.(72), 


